Kinetic Monte Carlo simulations of the growth of polymer crystals 
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Based upon kinetic Monte Carlo simulations of crystallization in a simple polymer model we 
present a new picture of the mechanism by which the thickness of lamellar polymer crystals is 
constrained to a value close to the minimum thermodynamically stable thickness, Zmin- The free 
energetic costs of the polymer extending beyond the edges of the previous crystalline layer and of a 
stem being shorter than /min provide upper and lower constraints on the length of stems in a new 
layer. Their combined effect is to cause the crystal thickness to converge dynamically to a value close 
to Zmin where growth with constant thickness then occurs. This description contrasts with those 
given by the two dominant theoretical approaches. However, at small supercoolings the rounding of 
the crystal profile does inhibit growth as suggested in Sadler and Gilmer's entropic barrier model. 
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I. INTRODUCTION 

On crystallization from solution and the melt many 
polymers form lamellae where the polymer chain tra- 
verses the thin dimension of the crystal many times fold- 
ing back on itself at each surface.1^1 Although lamellar 
crystals were first observed over forty years ago their 
physical origin is still controversial. It is agreed that the 
kinetics of crystallization are crucial since extended-chain 
crystals are thermodynamically more stable than lamel- 
lae. However, the explanations for the dependence of 
the lamellar thickness on temperature offered by the ta'Q 
dominant theoretical approaches appear irreconcilable E'cl 
The lamellar thickness is always slightly greater than 
^minj the minimum thickness for which the crystal is ther- 
modynamically more stable than the melt; Imin is approx- 
imately inversely proportional to the degree of supercool- 
ing. 

The first theory, which was formulated by Lauritzen 
and Hoffman (LH) soou, after the initial discovery of the 
chain- folded crystals ,LrO invokes surface nucleation of a 
new layer on the thin side faces of the lamellae as the 
key process. It assumes that there is an ensemble of 
crystals of different thickness, each of which grows with 
constant thickness. The crystals which grow most rapidly 
dominate this ensemble, and so the average value of the 
thickness in the ensemble, which is equated with the ob- 
served thickness, is slightly larger than the thickness for 
which the crystals have the maximum growth rate. The 
growth rates are derived by assuming that a new crys- 
talline layer grows by the deposition of a succession of 
stems (straight portions of the polymer that traverse the 
crystal once) along the growth face. The two main factors 
that determine the growth rate are the thermodynamic 
driving force and the free energy barrier to deposition of 
the first stem in a layer. The former only favours crys- 
tallization when the thickness is greater than Imin] the 
latter increases with the thickness of the crystal because 
of the free energetic cost of creating the two new lateral 



surfaces on either side of the stem and makes crystalliza- 
tion of thick crystals prohibitively slow. Therefore, the 
growth rate passes through a maximum at an intermedi- 
ate value of the thickness which is slightly greater than 

The second approach, which was developed by Sadler 
and Gilmer and has been termed the entropic barrier 
model, is based uawi the interpretation of kinetic [Maate 
Carlo simulationsQu and rate-theory calculationsETM of 
a simplified model of polymer crystal growth. As with 
the surface nucleation approach, the observed thickness 
is suggested to result from the competition between a 
driving force and a free energy barrier contribution to 
the growth rate. However, a different cause for the free 
energy barrier is postulated. As the polymer surface in 
the model can be rough, it is concluded that the details 
of surface nucleation of new layers are not important. In- 
stead, the outer layer of the crystal is found to be thinner 
than in the bulk; thisjrounded crystal profile prevents fur- 
ther crystallization.1121 Therefore, growth of a new layer 
can only begin once a fluctuation occurs to an entropi- 
cally unlikely configuration in which the crystal profile is 
'squared-off '. As this fluctuation becomes more unlikely 
with increasing crystal thickness, the entropic barrier to 
crystallization increases with thickness. 

Although both approaches are able to describe cor- 
rectly some of the basic phenomenology of polymer crys- 
tallization, both have questionable aspects. As Frank 
and Tosi pointed out, one implication of the LH assump- 
tions is that the thickness of an individual crystaLshould 
not vary even when the temperature is changed.c3 How- 
ever, this is clearly contradicted by experiment: when the 
temperature is changed the thickness of a growing crys- 
tal adjusts U^iAe new temperature producing a step on 
the lamella.ll3lij Sadler has also extensively questioned 
the adequacy of the LH approach, in particular to ex- 
plain curved crystal, habits and the effect of twins on 
the growth rate.t5llj Furthermore, there is the so-called 
'5Z-catastrophe': at too large supercoohngs the predicted 
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thickness goes to infinity unless a parameter in the model, 
the apportionment factor 4", is rather arbitrarily set to 
zero.B 

In the entropic barrier approach it is not clear whether 
approximations such as the implicit representation of the 
chain connectivity and chain folds by a set of simple 
growth rules and the neglect of the energetic contribu- 
tion of chain folds to the free energy of the fold surface 
compromise its conclusions. Furthermore, although the 
simulations have been interpreted in terms of an entropic 
barrier, no direct evidence of this free energy barrier has 
been provided.E3 

There have been a number of works which have taken 
the basic LH approach but relaxed some roi the con- 
straints. For example- both Frank and Tosilla and Lau- 
ritzen and PassagliaE^ have allowed the stem lengths 
within a layer and between layers to vary. Later, 
Point allowed the first stem to be deposited in steps 
rather than as a G«ijiiplete stem; this change prevents 
the (5Z-catastrophe.E3H In the same spirit DiMarzio and 
Guttman wanted to allow every stem to be deposited 
in steps and to allow a fold to be formed at any stepi 
however they could only consider very restricted cases. cil 
The approach used in all these papers was to attempt to 
find analytical solutions to the rate-theory equations, but 
as the number of possible pathways considered increased 
the problems became increasingly intractable. The nat- 
ural solution to this difficulty is to use computer simu- 
lation techniques such as kinetic Monte Carlo (KMC). 
However, all these papers are over fifteen years old, and 
at that time such an approach was not so feasible; only 
Point made an attempt to apply computational methods 
to the problem.Ea 

Therefore, a thorough application of computer simu- 
lation to an unrestricted version of the LH approach is 
long overdue. In this paper we do just this. We ex- 
amine the growth of new crystalline layers from solution 
on a surface that represents the growth face of a poly- 
mer crystal when there are no constraints on the stem 
length or on when a fold can form. Our model is, in 
some ways, midway between the original LH model and 
the model used by Sadler and Gilmer. In the latter, a 
realistic representation of the polymer is sacrificed in or- 
der to be able to consider, for example, rough growth 
surfaces. Consequently, our work, as well as providing 
a test of the two theories, could play an important role 
in understanding the relationship between them. This 
work should also be seen in the context of an increasing 
application of simulatioa, to, study processes relevant to 
polymer crystallization.t3~E3 



II. METHODS 

In our simulations the polymer is represented by a self- 
avoiding walk on a simple cubic lattice. There is an at- 
tractive energy, — e, between non-bonded polymer units 



on adjacent lattice sites and an energetic penalty, eg, for 
each kink (a 'gauche bond') in the chain. As the surface 
in our system represents the growth face of a polymer 
crystal the — e interaction also applies to polymer units 
in contact with the surface, e can be considered to be 
an effective interaction representing the combined effects 
of polymer-polymer, polymer-solvent and solvent-solvent 
interactions, and so we can use our model in a simpli- 
fied representation of the crystallization of a semi-flexible 
polymer from solution. The behaviour of the polymer is 
controlled by the ratio kT/e; large values can be consid- 
ered as either high temperature or good solvent condi- 
tions, and low values as low temperature or bad solvent 
conditions. The parameter eg defines the stiffness of the 
chain. The polymer chain is flexible at £^=0 and becomes 
stiffer as eg increases. Here, we use eg = 8e_however sim- 
ilar results are obtained at any positive egE3 

This interaction scheme has been recently used 
to investigate-, -tlie phase behaviour of isolated 
homopolymeia^jcj and a homopolymer in the presence 
of a surface.Lj The low temperature behaviour is of par- 
ticular interest since the polymers were found to adopt 
cuboidalE3 (rectangular when on a surfaceE3) 'crystalline' 
configurations involving chain folding. However, these 
structures were observed for thermodynamic reasons — 
they are lowest in energy — whereas chain-folding occurs 
in polymer crystallization for kinetic reasons. 

In our simulations we wish to investigate the growth 
of new crystalline layers on the growth face of a crystal. 
To achieve this we choose not to use conventional Monte 
Carlo or molecular dynamics simulations since it is not 
feasible with these methods to probe the desired tiroa 
and length scales. Instead we use kinetic Monte Carlo.Eil 
However, this, of course, means that we have to make 
certain assumptions about the processes contributing to 
polymer crystallization. In our model when a polymer 
crystallizes on the growth face it forms a two-dimensional 
crystalline layer. The crystals formed only have adjacent 
re-entry of the stems; i.e. only tight folds where the new 
stem is adjacent to the previous stem are allowed. An 
example configuration is shown in Figure |l]. 

To make the problem more computationally tractable 
we only model the crystalline portion of the polymer ex- 
plicitly. The rest of the chain is assumed to behave like an 
ideal coil. If we consider the coil to be a two-dimensional 
coil adsorbed onto the surface the free energy is 

^idcai,2D = -NcoiikT log (1 -h 2 exp{-f3eg)) - iVcoue, (1) 

whereas if the coil adopts a three-dimensional configura- 
tion in solution the free energy is 

^idoaUD - - A^coiifcT log (l-f 4 exp(-/3eg)) , (2) 

where A'coii is the number of units in the coil state. These 
expressions ignore any energetic contributions from con- 
tacts between polymer units in the coil. This approxi- 
mation becomes worse as the temperature decreases, be- 
cause the disordered chain would tbep be expected to 
have a dense collapsed conformation.l23 
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For a single homopolymer (with the current interac- 
tions scheme and with aU units exphcitly represented) 
in the presence of an infinite polymer-Hke surface, the 
lowest energy configuration is a two-dimensional crys- 
talline configuration. Simulations have shown that on 
melting of the crystalline polymer a two-dimensional dis- 
ordered state was formed for all values of eg, and only 
at higher temperature does the polymer become more 
three-dimensional.El This suggests that Equation (|l]) is 
an appropriate form to describe the coil. However, when 
considering growth on a thin lamella, the surface area 
on which the coil can absorb is much reduced. To stay 
on the surface the coil must adopt an anisotropic con- 
figuration in the plane which would reduce its entropy 
and thus reduce the free energetic advantage of a 2D coil 
compared to a 3D coil. Although we predominantly use 
Equation (|l|) to describe the coil, we also consider the 
effects of assuming the coil is 3D. 



FIG. 1. An example configuration for the crystalline por- 
tion of the polymer which shows the possible moves (rep- 
resented by arrows) in the kinetic Monte Carlo simulation. 
For this configuration, where growth is allowed to occur from 
either end, there are five possible moves. The dashed lines 
mark the boundaries of the underlying surface 

At each step changes of configuration can only occur 
at the ends of the crystalline portion of the polymer. The 
possible processes are for the crystalline part of the poly- 
mer to grow by one unit by the extension of one of the 
stems or by the formation of a gauche bond as part of a 
fold, or for the crystalline part of the polymer to shrink 
by one unit. To preserve adjacent re-entry after the for- 
mation of the first of the two gauche bonds of a fold the 
possible moves are just the completion of the fold and 
the removal of the first gauche bond. These possibilities 
are illustrated in Figure ||. Generally we allow growth to 
occur from both ends of the crystalline configuration, but 
occasionally, in order to compare with theories of poly- 
mer crystallization, we only allow growth to occur from 
one end. Physically, these two scenarios correspond to 
the initial crystal nucleus being in the mid-section of the 
polymer, or at the end of the polymer, respectively. We 



see no reason why the latter should always be the case. 
The choice only makes a difference to the growth rate, 
which is of course twice as fast when growth can occur 
from both ends, and the structure of the crystal near to 
the initial nucleus. 

Note that our scheme does not allow for the annealing 
of the part of the crystallite that has already been de- 
posited. This will not affect our conclusions concerning 
the average thickness of the deposited layers but it may 
lead to an unrealistically rough appearance of the crystal. 
We will come back to this issue later in the paper. 

The change in the free energy of the polymer on tak- 
ing a step, AA, is given by = Ai^xtai + AAcoii, where 
AiSxtai is the change in energy of the crystalline configu- 
ration and AAcoii is calculated from Equation (|I]) or (|^). 
A rate is assigned to each of the possible moves. Each 
move is assumed to be an activated process with a barrier 
that is an energy A above the higher free energy state. 
The prefactor is assumed to be the same for all processes, 
and we set it to 1 thus defining a time unit. Therefore, 
the rate for a step is given by 



fcy = exp ( ] min ( 1, exp 



kT 



(3) 



As the factor exp(— A/fcT) scales all the rates, simula- 
tions do not need to be carried out at different values of 
A, but the results for one value can be scaled onto any 
other. 

In the KMC simulation at each step we choose ran- 
domly one state, j, from those connected to the current 
state, i, with a probability given by 



p. 

J- 7,1 



'ij' 



and update the time by an increment 



At = - 



log(p) 

J2i ^ij 



(4) 



(5) 



where p is a random number in the range [0,1]. In this 
way the KMC algorithm simulates a stochastic process 
described by a Poisson distribution. 

As the expressions for the various rates do not depend 
on the length of the polymer we do not set a limit to 
the length of the polymer, but allow each new crystalline 
layer to be formed from a single polymer. 

Most of the results in this paper are for the growth 
of a single crystalline layer on a growth face which is 
of uniform thickness. However, we have also considered 
the situation where many layers are grown on the initial 
crystal. In these simulations all except the first layer are 
grown on top of the previously-grown layer and we only 
start to grow a new layer once the previous layer has 
been completed (when the crystalline polymer spans the 
periodic boundary conditions in the a-direction). 
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FIG. 2. Polymer configurations for a new crystalline layer grown at T = 2.75 ek~^ on a surface of thickness; (a) oo, (b) 
50 and (c) 21. Each configuration contains 150 stems and a dashed line marks the boundaries of the underlying surface; the 
3/-axis origin is at the centre of the surface. 




1000000 2000000 3000000 4000000 5000000 
steps 



The number of units in the crystalline 
as a function of the number of steps for (a) 



FIG. 3. 

layer, A''xtai, 

T = 3.7efc"\ (b) T = 4.1 £fc"\ The solid line (s) shows the 
behaviour of a single layer, whereas the dashed line (a) is 
an average over 50 layers. The underlying surface is infinite. 
Initially, there is in (a) a single unit on the surface and in 
(b) a crystal 100 stems wide, each 120 units long. 

The thermodynamics of this model can be easily 
worked out if we ignore any entropic contribution to the 
crystal from variations in the stem length. The energy of 
a crystal with constant thickness, /, is 



xtal 



2e„ 



2e 



(6) 



where N^tai is the number of units in the crystalline 
configuration. The positive term is due to the sur- 



face energy of the fold surface. The minimum stable 
thickness, /min, can be calculated through the equation 
E^tai/N^tui = Acoii/Ncoih This gives 



2e„ 



'min.2D — 



fcTlog(l + 2exp(-/3eg)) 



and 



2£g + e 

2e-fcTlog(l + 4exp(-/3eg))' 



(7) 



(8) 



The neglect of the entropy of the crystal means that these 
expressions are upper bounds. The temperature at which 
the surface of the polymer crystal loses its crystalline or- 
der, the 'melting point' T^, is given by the temperature 
for which the denominator is zero. 



III. RESULTS 

Some aspects of the typical behaviour of our model 
for the growth of a single layer on a surface are illus- 
trates in Figures ^ and ^. Below the melting point 
(r,„ = A.06 ek~^), for a sufficiently thick surface, a new 
crystalline layer grows. However, at low supercoolings 
the initial nucleation of a new layer can be quite slow. In 
the example shown in Figure there is no net growth 
in the first million time steps, and only once a viable 
nucleus has formed does growth occur with relatively lit- 
tle impediment. This example clearly shows the effect 
on the dynamics of the free energy barrier for the initial 
nucleation of a layer. The time scale for the initial nucle- 
ation event decreases with increasing supercooling and at 
large supercoolings it is not noticeable in the trajectories. 
(Very close to it is only feasible to grow new layers 
when a seed crystal is introduced.) 

Above T„i, as one would expect, if we start the sim- 
ulation with an initial seed crystal it duly dissolves. In 
the example shown in Figure ^b the temperature is just 
above Tm, and so the rate of dissolution is slow because 
of the nearly-flat free energy landscape. In every 1000 
steps there is only 3 more dissolution than growth steps. 
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In the LH theory the thickness of a layer is determined 
by the length of the first stem that is deposited. However, 
as the free energy barrier associated with the first stem is 
proportional to the length of the stem, one could imagine 
that growth would be more rapid for a layer where the 
stem length gradually increases to its average value as 
crystallization progresses than for a layer with constant 
stem length. In such a configuration the lateral surface 
free energy is paid for 'in installments' rather than all 
initially and growth is likely to be more rapid because of 
the lower initial free energy barrier. 
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growth which causes the stem length to be constant. (Al- 
though subsequent annealing of the new crystalline layer 
may reduce the roughness.) Secondly, we have tested 
whether there is any systematic dependence of the aver- 
age stem length on the position of the stem with respect 
to the initial nucleation site (Figure ^). Allowing the 
crystal to grow from only one end most closely mimics 
the conditions of the LH theory in which stems are grown 
one at a time. In this case we see that the suggestion we 
made above is confirmed. The length of the first stem is 
significantly shorter than the average stem length in the 
layer, and the stem length converges to the average value 
for the layer as the crystal grows. 

Interestingly when we allow growth to occur from both 
ends the initial behaviour is very different (Figure ^). 
This is because it is now possible to form an initial nu- 
cleus which consists of two stems connected by a fold 
which grow simultaneously, a scenario not considered in 
the LH theory. Such a nucleus is energetically favoured 
because of the interactions between the two stems. As 
it avoids the large fee energy barrier associated with a 
single-stem nucleug23 there is now no advantage in the 
initial stems being shorter than the average in the layer 
(Figure ^). The possibility of a two-stem nucleus was 
first suggested by Pointed and recently strong evidence 
for this nucleus has been obtained in a study of the fsee 
energy profiles along specific crystallization pathwaysE3 
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FIG. 4. The dependence of the average stem length on the 
distance of the stem from (a) the initial nucleation site, and 
(b) the centre of an initial crystal seed at T = 2.75 ek~^ for 
growth on an infinite surface. In (a) growth starts with a 
single polymer unit on the surface and we consider the cases 
where growth is allowed at one end or both ends of the crys- 
talline portion of the polymer. In (b) the crystal seeds are 
3 stems wide; the lengths of the stems in the seeds are as 
labelled. 

Our model allows us to examine these questions easily. 
Firstly, inspection of example configurations shows that 
there is considerable variation in the stem lengths within 
a layer (Figure ||). There is nothing in the kinetics of 



FIG. 5. The dependence of the stem length on the length 
of the previous stem for the growth of a single layer on a sur- 
face 50 units thick at T = 2.75 efc~^. 

This finding about the structure of the initial nucleus 
significantly dents the LH theory because its assumption 
about the nature of the initial nucleus, and the conse- 
quent free energy barrier, is crucial to the theory. How- 
ever, this finding is rendered somewhat irrelevant because 
the above results also undermine a more fundamental 
tenet of the LH theory, namely the assumption that the 
thickness of a layer is determined by the initial nucleus. 
Further confirmation that this assumption does not hold 
comes when we examine the growth from initial seed crys- 
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tals. Whatever the thickness of the initial seed the thick- 
ness of the growing crystal converges to the same value 
(Figure These results imply that the thickness of a 
crystalline layer must be determined by factors which are 
operating on the deposition of each stem and not those 
specific to the initial stems. We shall examine what these 
factors are in detail later. 

In Figure || we show the correlations between succes- 
sively grown stems. The figure has the form of a fixed- 
point attractor. It shows there is one stem length, I*, for 
which the average length of the next stem is the same as 
the previous (for Figure || r=36), and this corresponds 
to the average stem length in the layer. For stems longer 
(shorter) than I* the length of the next stem is on aver- 
age shorter (longer) bringing the stem length nearer to 
I*. Frank and Tosi have previously drawn a similar con- 
clusion from their extension of the LH approach in which 
a single change in the stepi length is allowed during the 
deposition of a new layer .E3 
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FIG. 6. (a) The average stem length and (b) the growth 
rate (polymer units per exp( — A/fcT) units of time) in the 
new crystalline layer as a function of temperature. The dif- 
ferent curves are for different thicknesses of the underlying 
surface, as labelled. 



25 30 35 

.surface thickness 




35 40 
surface thickness 

FIG. 7. The average stem length in the new crystalline 
layer as a function of the thickness of the underlying surface 
for (a) r = 2.75efc"\ (b) T = 3.375 efc~\ The dotted lines 
show how the thickness changes on addition of successive 
layers to a 50-unit thick surface. 



In Figure |^a we show how the thickness of a new layer 
depends on the temperature and the thickness of the 
growth surface. It can be immediately seen from the 
fact that the lines are not horizontal that the thickness 
of a new layer is not necessarily the same as the previous 
layer, thus contradicting another of the assumptions of 
the LH theory. Instead the temperature dependence of 
the curves reflects the factors that determine the thick- 
ness of a new layer. All the thickness curves rise as the 
temperature approaches Tm because of the rise of Zmin; 
they end at the temperature where it is no longer pos- 
sible to grow a new layer which is thermodynamically 
stable. (Figure ^ shows that the growth rate of the 
crystal goes to zero at these end points.) At low tem- 
perature the thickness also increases, in this instance, 
because it becomes increasingly |difficult to scale the free 
energy barrier for forming a foldE3 and so on average the 
stem continues to grow for longer. However, except for 
growth on an infinite surface, this rise is checked by the 
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presence of the edge of the underlying surface. It is un- 
favourable for the polymer to overhang the edge because 
these units do not gain the — e energy of interaction with 
the surface. 

An understanding of how the thickness of a polymer 
crystal is determined can be gained from these results but 
to do so it is better to plot the thickness of the new layer 
against the thickness of the growth surface at constant 
temperature, as in Figure |^. In our KMC simulations we 
do not allow the crystalline layer to anneal. However, it 
is likely that some annealing does take place before the 
next layer is deposited. In what follows we consider two 
limits. The first is complete annealing. In that case we 
may replace the previous crystalline layer by a flat sub- 
strate with the same thickness. The second limit is no 
annealing. This we discuss with reference to Figure 
Let us first focus on the complete annealing limit. In 
that case we can use Figure ^ to consider the effects of 
growing successive layers on top of each other. Follow- 
ing the dotted lines shows one what would happen for 
growth on an initial surface 50 units thick. For example, 
at T = 2.75ek~^ the first layer is 36 units thick, the sec- 
ond 28, the third 23, . . .The thickness converges to the 
value I** at which the curve crosses y = x (/**=21 at 
T ~ 2.75 e/c^^), i.e. to the point where the thickness of 
the new layer is the same as the previous, and then the 
crystal continues to grow at that thickness. Figure ^ is 
another example of a fixed-point attractor. 

Again this picture is very different from that presented 
by the LH theory. It shows that it is inappropriate to 



compare the growth rates of crystals of different thick- 
ness because there is only one dynamically-stable thick- 
ness, I**, for which growth at constant thickness can oc- 
cur. For crystals which initially have a thickness differ- 
ent from I**, during growth the thickness. will converge 
to I**, as has been observed in experiment .EjO It is also 
interesting to note that the growth rate at I** is not a 
maximum . Quite on the contrary, the growth of a thick 
crystal slows as the thickness decreases to I** (Figure ||). 
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FIG. 8. The growth rate (polymer units per exp(— A/fcT) 
units of time) of the new layer as a function of the thickness 
of the underlying surface for T — 2.75 ek~^ . A dotted verti- 
cal line has been placed at the thickness, l" . 
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Although this viewpoint contrasts with the LH theory, 
it was realizedr4B|--a|Winiber of the early polymer crystal- 
lization paperglaO'Ea that stable growth could only occur 
at the thickness for which a new layer has the same thick- 
ness as the previous. However, since then this insight has 
for the most part been neglected, and, to the best of our 
knowledge, the iterative maps that underlie it have not 
been previously visualized. Also, in recent simulations by 
Chen and Higgs the lamellar thickness was observed to 
converge to the same value-|irrespective of the thickness 
of the initial seed crystal,E3 but the significance of this 
behaviour was not fully realized. 

We can better understand the reasons for this be- 
haviour by examining representative polymer configura- 
tions (Figure |^) and the probability distributions of the 
stem length (Figure ||) for the growth of a single layer 
at a number of surface thicknesses, /min places one con- 
straint on the stem length; only a small fraction of the 
stems can be shorter than if the layer is to be ther- 
modynamically stable. The boundary of the growth face 
places the second constraint on the stem length; it is 
energetically unfavourable for the polymer to extend be- 
yond the edges of the underlying surface. Even in the 
absence of this constraint, i.e. on the infinite surface, the 
stem length remains finite because at every step theceis. 
always a finite probability that a fold will be formedjH3'0 
and so the probability distribution decays to zero at large 
stem lengths (Figure ||a). For the infinite surface any 
stem that is longer than ^min is stable with respect to 
the coil and so a very wide range of stem lengths is ob- 
served (the standard deviation in the stem length is '^35 
at r = 2.75efc-i). 

As can be clearly seen in Figures |^ and ||, the thickness 
of the growth face exerts an increasingly strong influence 
on the new layer as the thickness decreases. For a sur- 
face 100-units thick the probability distribution of the 
stem length is very similar to the infinite case — again 
the maximum occurs at '^40-50 — except that the long 
tail of the distribution is cut off at 100 because the poly- 
mer rarely extends beyond the edge of the growth face. 
At this thickness the probability of the stem length be- 
ing greater than the surface thickness is much less than it 
being smaller and therefore, the new layer is significantly 
thinner than the underlying surface. As the surface thick- 
ness decreases the probability distributions of the stem 
length becomes increasingly narrow and the difference in 
probability between the stem length being greater or less 
than the surface thickness diminishes. This is because 
the range of viable stems — those with lengths larger than 
^min but which do not significantly overhang the edge of 
the surface — become smaller. Finally, at I** as the sur- 
face thickness approaches ^min the probability distribu- 
tion become symmetrical about the surface thickness and 
the thickness of the new layer becomes equal to the thick- 
ness of the growth surface (Figure |^). For these reasons 
the thickness at which stable growth occurs is close to 
^min- For thicknesses less than /** the asymmetry of the 
probability distribution is reversed (Figure |9|f ) . 



0.06 




-60 -40 -20 20 40 60 



FIG. 10. The probability distribution for tfie position of 
the ends of the crystalline portion of the polymer for growth 
of a single layer on surfaces of different thickness, as labelled. 
The centre of the surface is at 0, and dotted lines have been 
included at the positions of the boundaries of the surface. 
T = 2.75 efc"\ 



It is also possible to understand why the dependence 
of the growth rate on thickness (Figure ^) differs from 
the predictions of the LH theory. The cause is not the 
decrease of the thermodynamic driving force as the thick- 
ness decreases (this causes the growth rate to go to zero 
at Imin) since it is common to both this model and the 
LH theory. Rather, the difference arises mainly from the 
LH assumption that complete stems are deposited in a 
single step. Based upon comparisons of this approach 
with a fine-grained kinetic description of the deposition 
of a stem, Frank and Tosi showed that the LH complete- 
stem approach is reasonable if the only point on the 
fine-grained pathway at which a fold can occur is once 
a stem is complete. However, as has been pointed out 
previously, this condition cLqcs not hold for more realistic 
multi-pathway approaches J3 That the polymer can form 
a fold at any point and can overhang the edge of the crys- 
tal in our model has a considerable effect on the growth 
rates. Such processes can lead to configurations that are 
not therniodynamically viable, such as when a new stem 
starts to form even though the previous is less than /„ii,i 
or when a stem significantly overhangs the edge of the 
crystal. These configurations have to be removed before 
growth can continue. The effect of such 'blind alleys' is 
to reduce the growth rate. As the thickness of the surface 
decreases the range of stem lengths which give viable con- 
figurations decreases and so an increasing amount of time 
is spent searching blind alleys. These effects are reflected 
in the probability distributions for the position of the end 
of the crystalline polymer (Figure |l0|). As the thickness 
decreases the end spends an increasing proportion of the 
time near the edge of the growth surface waiting for the 
polymer to fold at a point when the stem has a length 
that is viable for further growth. 

At T = 2.75 ek^^ there is no barrier to growth due to a 
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rounded crystal profile. When the thickness is I** the new 
layer has the same thickness as the previous layer. How- 
ever, this scenario does not hold for all temperatures. At 
any temperature there is a thickness below which a new 
layer cannot grow because the surface is too thin, (e.g. 
the last point of the curve in Figure 0a is for a surface 
19 units thick) and there is no a priori reason why this 
must occur after the thickness curve has crossed y = x. 
Indeed, for T > 3.2efc~^ there is no thickness for which 
successive layers have the same thickness. For example, 
aX T — 3.375 efc~^ after the growth of two layers on a 
50-unit thick surface the outer layer is ~29 units thick 
(Figure the crystal then stops growing because the 
outer layer is too thin. For these smaller supercoolings, 
as suggested in the entropic barrier model, the rounding 
of the crystal profile inhibits growth. 

To overcome this barrier requires a cooperative mech- 
anism whereby a new layer takes advantage of (and then 
locks in) dynamic fluctuations in the outer layer to larger 
thickness. The presence of such fluctuations is shown in 
the probability distributions of Figure]^; however, growth 
stops in our model because we attempt to grow a new 
layer on an outer layer that is static. As overcoming the 
barrier would be most rapid when the magnitude of the 
fluctuations is the minimum necessary, we expect that 
this mechanism would lead to the crystal continuing to 
grow with the smallest thickness for which a new layer 
can grow, e.g. at T = 3.375 ek~^ this is a thickness of 30 
(Figure ^). This mechanism again leads to a thickness 
for the polymer crystal which is close to /min- 




FIG. 11. Cut through a polymer crystal which was pro- 
duced by the growth of twenty successive layers on a surface 
with a uniform thickness of 50 units at T = 2.0efc~^. The 
stems are represented by vertical cuboids. The cut is 16 
stems wide. 

It is interesting to note that behaviour similar to our 
results has been found in some of the extensions of the 
LH theory which include multiple pathways in some re- 



stricted or approximate manner. In a typically prescient 
paper Frank aed Tosi found that growth ceases at low 
supercoolings;ll3 however, without a mechanism to allow 
growth to continue, they concluded that polymer crys- 
tallization experiments could not correspond to this tem- 
perature range. Lauritzen and Passaglia were also aware 
that growth would cease in their model at low supercool- 
ings, so they introduced an ad hoc eneMetic term into 
their rate constants that prevented this.til'Ej 

In using Figure |^ to consider the growth of successive 
layers on the initial growth face we are assuming that all 
the roughness in the position of the fold surface is an- 
nealed out before a new layer starts to form. To check 
that the behaviour we observe is not dependent on this 
assumption we have also performed simulations where we 
grow new layers directly on top of the previous grown lay- 
ers. A typical cut through a crystal that was produced 
by this method is depicted in Figure |ll|. We see the same 
mechanism of thickness determination as before: within 
5-10 layers the thickness of the crystal converges to I** 
and then growth continues at that thickness. 

H 



j j+l 



FIG. 12. An example of a mechanism by which the 
roughness in the crystalline layer could be annealed. The 
stem j slides upwards decreasing the length of stem j — 1 
by one and increasing the length of stem ji -I- 1 by one. The 
change in energy for this rearrangement is — e. 

It is not clear which of these two limiting approaches 
is more realistic. In an atomic-force microscopy (AFM) 
study of the fold surface of polyethylene only on «, m'-. 
nority of crystals could the folds be resolved clearlyEJ'EEl 
suggesting that for most crystals the fold surface retains 
some roughness. However, the crystalline configurations 
generated by our model when the thickness of the growth 
face is much larger than Imin (Figure |^) are clearly too 
rough. One would imagine that local rearrangements 
such as the one depicted in Figure |lj would act to reduce 
the roughness whilst keeping the average stem length in 
the layer constant. More uniform configurations could 
have been obtained by including these mechanisms in 
our model but we did not do so because we would have 
to have set the relative rates of growth and annealing to 
an arbitrary value. To study the multi-layer growth ver- 
sion of our model more systematically we would probably 
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need to include these mechanisms because the roughness 
can at times hinder the growth; e.g. it can be hard for a 
new layer to grow across a section of the previous layer 
that involves a particularly short stem. Chen and Higgs 
experienced similar difficulties in their simulations.EZI 
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FIG. 13. The dependence of the average stem length in 
the new crystalline layer on the temperature and thickness 
of the surface when the disordered polymer is assumed to 
be an ideal three-dimensional coil. In (a) the labels give the 
value of the thickness of the underlying surface, and in (b) 
the value of the temperature in ek~^ . 

To determine the effect of our assumption that the 
disordered polymer is adsorbed onto the surface we have 
also performed simulations where the coil is assumed to 
be three-dimensional. Figure [l^ a shows the dependence 
of the thickness of a new layer on the temperature and 
thickness of the surface. The figure has a similar form 
to Figure ^a except that the rise in the layer thickness 
at low temperature is not checked by the thickness of 
the underlying layer. The free energy difference between 
the coil and the crystal at low temperature is now large 
enough to allow the polymer to extend beyond the edge 
of the growth face. This has a corresponding effect on 
I** . When the coil was assumed to be two-dimensional 



I** was a monotonically increasing function of tempera- 
ture which was always slightly larger than /min- However, 
now the value of I** increases at low temperatures; e.g. at 
T = 3.0efc-i r* = 21 whereas at T = 3.5 efc-^ 1** = 18 
(Figure [T^b). Assuming the coil to be three-dimensional 
does not prevent the rounding of the crystal profile which 
causes growth to stop in our model at low supercoolings; 
however, the range of temperature for which this scenario 
holds is reduced (T > 3.8eA:~^) because of the greater 
thermodynamic driving force for crystallization. 



IV. CONCLUSIONS 

In summary, our results present a new picture of the 
mechanisms that cause the thickness of lamellar polymer 
crystals to be constrained to a value just above Zmin. The 
free energetic costs of the polymer extending beyond the 
edges of the previous crystalline layer and of a stem being 
shorter than Zmin provide upper and lower constraints on 
the length of stems in a new layer. Their combined effect 
is to cause the crystal thickness to converge to a value 
close to Zmin whcrc growth with constant thickness then 
occurs. 

This convergence of the thickness has been observed 
in experiments in which the thickness of growing poly- 
mer crystals adjusts to a change in temperaturell3'Efl and 
in which lamellae form by epitapdal crystallization onto 
extended-chain crystalline fibresEj. It would be very in- 
teresting if AFM could be used to probe the profiles of 
the steps on the lamellae that result from temperature 
changes. From these profiles it would be possible to con- 
struct maps similar to Fig. ^ which could confirm the 
crucial role of the thickness of the growth face on the 
properties of a new crystalline layer. AFM could be also 
used to study the profile of the crystal close to the growth 
face to examine whether rounding of the crystal edge oc- 
curs. 

Our results have significant implications for the LH 
surface nucleation theory. Although in some ways our 
model is in the LH tradition — new layers nucleate and 
grow by the addition of a succession of stems along the 
growth face — the removal of many of the LH constraints 
leads to a significantly different picture for the mecha- 
nism that determines the thickness of polymer crystals. 
Our work undermines many of the assumptions that are 
crucial to the LH theory: the initial nucleus is not a single 
stem; the initial nucleus does not determine the thickness 
within a layer; the thickness of a new layer is not neces- 
sarily the same as the previous layer; the observed thick- 
ness corresponds to the one value of the thickness, /**, for 
which crystals can grow with constant thickness and not 
to the thickness of those crystals which grow most rapidly 
in a fictitious ensemble of crystals of different thickness 
that grow with constant thickness; the growth rates are 
significantly affected by 'blind alleys' (not included in the 
LH theory) which lead to configurations that are not vi- 
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able for further growth. Of course, the attractive feature 
of the LH theory is that it provides analytical equations 
for the thickness and growth rate of polymer crystals that 
are in reasonable agreement with experiment. However, 
our results clearly indicate that the LH equations should 
be viewed as no more than phcnomcnological because 
their derivation relies upon an incorrect description of 
the physics of polymer crystallization. 

The picture of polymer crystallization that comes from 
our results has more in common with the entropic barrier 
approach, and these two models are, in some ways, com- 
plementary. In the simulations of Sadler and Gilmer a 
realistic description of the polymer connectivity and folds 
was sacrificed in order to be able to study the full three- 
dimensionality of the problem, whereas here we have a 
more realistic description of the polymer but could not at 
the same time include the cooperative interlayer dynam- 
ics that are necessary for growth at low supercoolings. 
Although Sadler and Gilmer's simulations have been in- 
terpreted in terms of the effects of a competition between 
an entropic barrier and a thermodynamic driving force 
on the growth rate, we expect that the mechanism of 
thickness determination is in fact similar to that for our 
model, namely that during gcowth the crystal thickness 
converges to a thickness The cessation of growth 

that we see in our model at low supercoolings confirms 
that the rounding of the crystal profile does play a role in 
polymer crystallization. However, in our model this ef- 
fect only occurs for a limited temperature range, whereas 
in the Sadler-Gilmer model rounding of the crystal profile 
occurs for all supercoolings. 

We believe that the results in this paper provide new 
insights into polymer crystallization, especially into the 
mechanism by which the thickness of a polymer crystal 
is attained. However, we stress that the model we use is 
simplified in order to make it computationally tractable. 
For example, we assume that nucleation only occurs once 
within a particular polymer; however, some simulations 
of a single polymer have seen crystalline domains develop 
within different parts of the polyme r wh ich at a later 
time coalesce to form a single crystal.ESES Furthermore, 
we only allow growth processes to occur within the plane 
of the growth face and only within the outer layer; as 
explained earlier, the cessation of growth at low super- 
coolings is a result of these simplifications. This study, 
then, represents only a starting point in the application 
of computer simulation to understand polymer crystal- 
lization. It is hoped that simulation will lead to an in- 
creasingly refined picture of the microscopic mechanisms 
involved in the growth of polymer crystals. 
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